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A short review on charmed and beauty hadron production in the lepton deep inelastic 
scattering off proton, in proton-proton and proton-antiproton coUisions at high energies is 
presented. It is shown that the existing theoretical and experimental information on charmed 
and beauty quark distributions in a proton is not satisfactory. A some procedure to study 
these distributions at LHC energies is suggested. 
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I. INTRODUCTION 

Various approaches of the perturbative QCD including the next-to-leading order calculations 
have been applied to construct the distributions of quarks in a proton. The theoretical analysis of 
the lepton deep inelastic scattering (DIS) off protons and nuclei provide rather realistic information 
on distribution of light quarks like u, d, s in proton. However, to find a believable distribution of 
heavy quarks like c(c) and especially b{b) in proton describing the experimental data on the DIS is 
non-trivial task. It is due mainly to small values of D- and B- meson yields produced in the DIS 
at existing energies. Even at the Tevatron energies the B- meson yield is not so large. However, 
at LHC energies the multiplicity of D- and B- mesons produced in p — p collisions will be larger 
significantly. Therefore one can try to extract a new information on the distribution of these heavy 
mesons in proton. In this paper we suggest some proposal to study the distribution of heavy quarks 
like c(c) and b(b) in proton based on the analysis of the future LHC experimental data. 
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II. HARD PARTON SCATTERING MODEL 

Usually, the multiple hadron production in hadron-nucleon collisions at high initial energies and 
large transfers is analyzed within the hard parton scattering model (HPSM) suggested in Refs. 

m 

and RefQ]. For example, the inclusive spectrum of hadron produced in the hard p — p interaction 
is presented in the following convolution form: 

where Eh,ph are the energy and three-momentum of the produced hadron h; Gqj^g{xi, kit) is the 
distribution of quark with flavor / or gluon in the first colliding proton depending on the Feynman 
variable xi and the transverse momentum pu; Gqj-^g{x2,k2t) is the same quark or gluon distribu- 
tion in the second colliding proton; s, t and u are the Mandelstam variables corresponding to the 



colliding quarks; da{s,i)/di is the differential 
z is the fraction of the produced hadron h 



cross section of the elastic parton-parton scattering; 



HQ. 

The HPSM has been improved significantly applying the QCD parton approach implemented in 
the modified minimal-subtraction renormalization and factorization scheme. It provides a rigor- 
ous theoretical framework for a global data analysis. In this framework there are two distinct 
approaches for next-to-leading order (NLO) calculations in perturbative QCD. 



The first calculation scheme is the so-ca^ 



led massive scheme or fixed-flavor-number scheme 



(FFNS) developed in Refs. and Refs..[6|, lZI]. In this approach the number of active flavors in 
the initial state is limited to nj = 4, e.g., u{u),d{d), s{s) and c(c) quarks being the initial partons, 
whereas the b(b) quark appears only in the final state. In this case the beauty quark is always 
treated as a heavy particle, not as a parton. In this scheme the mass of heavy quarks acts as cutoff 
parameter for the initial- and final-state collinear singularities and sets the scale for perturbative 
calculations. Actually, the FFNS with Uf = 4 is limited to a rather small range of transverse 
momenta pt of produced D or i?-mesons less than the masses of c or 6 quarks. In this scheme the 
m'^^/Pf terms are fully included. 

The another approach is the so-called zero-mass variable-flavor-number scheme (ZM-VFNS), 

riQ □ 

see Refs.[8|, lil] and Ref.|lCII] and references therein. It is the conventional parton model approach, 
the zero-mass parton approximation is applied also to the b quark, although its mass is certainly 
much larger than the asymptotic scale parameter Aqcd- In this approach the b{b) quark is treated 
as incoming parton originating from colliding protons or proton-antiproton. This approach can 
be used in the region of large transverse momenta of produced charmed or beauty mesons, e.g., 
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at pt > rric^h. Within this scheme the terms of order rri^jj/Pt can be neglected. Recently the 
experimental inclusive pt- spectra of i?-mesons in p—p collisions obtained by the CDF Collaboration 
[13I at the Tevatron energy ^/s=1.96 TeV in the rapidity region — 1 < y < 1 have been described 



ll| at > 10{GeV/c) using the non- 



rather satisfactory within this ZM-VFNS approach in Ref. 
perturbative structure functions. In another kinematic region, e.g., at 2.5{GeV/c < pt ^ 10{GeV/c) 
the FFNS model mentioned above allowed to describe the CDF data without using fragmentation 
functions of b quarks to i?-mesons. Both these schemes have some uncertainties related to the 
renormalization parameters. 

Looking the CDF data on the pt-spectra of i?-mesons produced in p—p collision at the Tevatron 

12 ] and Ref. [l^ one can find a difference between them especially at high 



energies published in Ref. 



values of pt- Therefore it would be very useful to have more precise data at different LHC energies. 

The experimental inclusive pt-spectra of D-mesons obtained by the CDF Collaboration at the 
same Tevatron energy and rapidity region [14] at pt > 5 (GeV/c) were described satisfactory within 
the so-called general-mass variable-flavor-number scheme (GM-VFN) [iS^ not assuming the zero 
mass for c quarks. 

To calculate the inclusive spectrum of hadrons, for example heavy mesons exploring eq.([T]) 
we have to know the distributions of quarks and gluons and their fragmentation functions (FF). 
Usually, they are calculated within the QCD using the experimental information from the DIS of 
leptons off protons and in the e"*" — e~ annihilation. The information on the gluon distribution 
and its fragmentation function can_be extracted from the experimental data on the jet production 
[1^ and their theoretical analysis 



13, 



18l |. Unfortunately the theoretical QCD calculation of the 
jet production has different sources of uncertainty. As is shown in Ref.[l(^] the main contribu- 
tion comes from uncertainty on parton distribution functions (PDFs) and is computed within the 
method suggested in Ref.|l9|. At low transverse momenta of jets pj.^* the uncertainty is small and 
approximately independent on the jet rapidity y-^^*. The uncertainty increases as p^^* and | y-'^* | 



increase. It can become about 130 percents 



Iq . To analyze the jet- and heavy quark-production at 



low Pt and large y or large values of Feynman variable x p one can apply the another nonperturbing 
QCD model, the so-called Quark-Gluon String Model. (QGSM). 



III. THE QUARK-GLUON STRING MODEL 



The QGSM is based on the 1/A^ expansion in QCD suggested in Refs.[20|, l2l|] instead of Og 
expansion that has the infrared divergence problem at 0, here is the number of flavors 
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or colors. The relation of the topological expansion over 1/N of the hadron-hadron scattering 



22 



23]. This 



approach 



amplitude to its f-channel one over Regge poles has been suggested in Refs. 
has been applied to a naly ze soft hadron processes at high energies, see for example Ref. 

It has been shown [24] that the main contribution to the inclusive spectrum of hadron produced 
m. p — p collisions at high energies comes from the so-called cylinder graphs corresponding to 
one-Pomeron and multi-Pomeron exchanges which are presented in Fig.l. 
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Fig. 1: The cylinder graph corresponding to the one-Pomeron exchange (left diagram) and the multi-cylinder 
graph corresponding to the multi-Pomcron exchanges (right diagram) . 

According to the QGSM, between quark q (diquark qq) and diquark qq (quark q) for colliding 
protons (Fig.l, left graph) the colorless strings are formed, then after their brake qq pairs are 
created which fragmentate to the hadron h. The right diagram in Fig.l corresponds to the creation 
of two colorless quark-diquark and diquark-quark strings and 2{n — 1) chains between see quark 
(antiquark) and antiquark (quark). The inclusive spectrum p/^ of hadron h produced \up—p collision 



corresponding to the one-Pomeron graph (Fig.l, left graph) has the following form 25l. l26l|: 



Ph{x,Pt) = aiF^{x+,pt)F^^{x.,vt)/F^^{{),pt) + F^^{x+,pt)F^{x^,pt) / F^{Q,pt) , (2) 

where o"i is the cross section of the 2-chains production, corresponding to the s-channel discon- 
tinuity of the cylinder (one-Pomeron) graph, usually it is calculated within the quasi-eikonal ap- 



proximation 



Fq[qq){x±,Vt) = ^ / ^a^l / ^ C^^P2t (a^l , Pit ) G^(gg) ( — , P2t )'^^^^ (Plt + P2t " Pt ) , (3) 

p, J x+ J Xl 

flavors ^ 

Here fq(gq){xi, pit) is the quark or diquark distribution function depending on the Feynman variable 
Xl and the transverse momentum of quark or diquark, z^^G'^^^^s^{z,p2t) = D^^^^s^{z,p2t) is the 
fragmentation function of quark or diquark to the hadron h; x± = ^{\/ Xf + x"^ it x) and xt = 
2m/it/-v/s, nifit = \/iT^'h + Pt^ ^ is the square of the initial energy in the p — p c.m.s. Actually, the 

gigq) ^ 

1, 



interaction function F^f^^^-^{xj^,pt) corresponds to the fragmentation of the upper quark (diquark) 
to the hadron /i, whereas F^ ,Jx-,pt) corresponds to the fragmentation of the down diquark 
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(quark) to h, see Fig.l (left diagram). The expression for the inclusive spectrum of the hadron h 
produced va. p — p collision corresponding to the right graph of Fig. 2 has more complicated form, 
see the details in Refs. 



25 



261]. 



All the quark distributions and fragmentation functions are related to the intercepts and slopes 
of Regge trajectories. For example, the distribution of c(c) quarks in a proton obtained within the 
QGSM and the Regge theory has the following form [25i] : 

An) ^ ,_,rc-"V.(0)fl _ ^a^^(0)-2aJv(0)+(ai^(0)-av,(0)+n-l 

where n is the number of Pomeron exchanges, aij(O), aAf(O) are the intercepts of the Reggeon and 
nucleon Regge trajectories, a^(0) is the intercept of the t/'-Regge trajectory, (5c(c) is a probability 
fraction of c(c) pairs in a quark see of the proton. 

The intercepts a_R(0) and aAf(O) are known very well from the experimental data on the soft 
hadron processes and the corresponding Regge trajectories have linear behavior as a function of 
the transfer t. As for the Regge trajectory a^(t), the information on its t-dependence is rather 
poor. As a function of t, it can be linear or nonlinear. 

The h{b) quark distributions and the fragmentation functions of h{h) to S-mesons can be also 
obtained within the Regge theory and the QGSM. They are related to the conventional Regge 
trajectories of light mesons and the Regge trajectory of the T-meson consisting of hh pair. The 



information on the t-dependence of the T-Regge trajectory is a.. 



The modified version of the QGSM suggested in Refs. 



251, 



so uncertain. 



261 ] including the transverse momen- 



tum dependence of the interaction functions Fqqg{x±,pt) allowed to describe the experimental data 
on inclusive spectra of L)-mesons produced in p — p collisions at the ISR energies and at moderate 
values of transverse momenta until ~ 4 — 5{GeV/c). However, the results are very sensitive to 
the value of the intercept a^(0). Note, that the QGSM in contrast to the perturbative QCD has 
not uncertainty related to the mass parameter. 



IV. PROPOSAL 



Concluding a short review on the HPMS and the QGSM we would like to propose to do a 
complex theoretical analysis both the jet production in p — p and p — p collisions within the 
perturbative QCD and the semi-hard production of D- and S-mesons in these reactions within the 
modified version of the QGSM. The semi-hard hadron processes mean the inclusive production of 
heavy mesons at not large transverse momenta pt and not small values of xp. We are going to 
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extract the gluon distribution in a proton and its fragmentation functions to heavy mesons from 
the QCD analysis of the jet production in p — p and p — p inelastic processes. Constructing a new 



modified version of the QGSM based on Refs. 



25 



we intent to include also the gluon-gluon and 



gluon-quark interactions using the obtained gluon distributions in proton and its FF. Then we will 
apply the suggested approach to describe all the existing experimental data on D- and S-mesons 
produced in p — p and p — p collisions at high energies and perform corresponding predictions at 
LHC energies. From this analysis we propose to extract a new information on the distribution of 
charmed and beauty quarks in proton. 
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